Pion production in the MiniBooNE experiment 
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Background: Charged current pion production gives information on the axial formfactors of nucleon resonances. 

It also introduces a noticeable background to quasi-elastic measurements on nuclear targets. 

Purpose: Understand pion production in neutrino interactions with nucleons and the reaction mechanism in 

nuclei. 

Method: The Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model is used for an investigation of neutrino- 
nucleus reactions. 

Results: Theoretical results for integrated and differential cross sections for the MiniBooNE neutrino flux are 
compared to the data. Two sets of pion production data on elementary targets are used to obtain limits for the 
neutrino-nucleus reactions. 

Conclusions: The MiniBooNE pion production data are consistent with the BNL elementary data if a small 
flux renormalization is performed while the ANL input data lead to significantly too low cross sections. A final 
determination of in-medium effects requires new data on elementary (p,D) targets. 
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I. INTRODUCTION 

In 2009 an unexpected experimental result appeared 
in studying quasielastic neutrino scattering from nuclei. 
The cross sections reported by the MiniBooNE collabo- 
ration at E v < 2 GeV were about 30% higher than 
the cross section measured by old bubble-chamber exper- 
iments in the 70s-80s. This excess was described in terms 
of a value for the axial mass that was significantly higher 
than the world-average value of about 1 GeV. At first 
sight, the new experiments have the advantage of huge 
statistics with millions of events recorded. One compli- 
cation, however, arises from the fact that these experi- 
ments all use nuclei as targets. All the measurements are 
thus influenced by nuclear effects. Indeed, the surprising 
result for the axial mass has since then found an expla- 
nation: The cross section labeled by the experiment as 
'quasielastic' contains in reality a significant admixture of 
2p-2h excitations that account for the difference between 
the data and theory using the standard value (ss 1 GeV) 
for the axial mass j. 

The quasielastic data themselves have not directly 
been measured but have been deduced from so-called 



quasielastic-like data by subtracting a background of 
'stuck-pion' events, i.e., events in which pions are first 
produced, but then reabsorbed again. This background 
was determined from calculations with an event genera- 
tor. Thus, the final QE + 2p-2h data invariably contain 
some model dependence. 

A broad comparison of generator predictions for QE 
and single-pion cross sections has shown major discrep- 
ancies (« 30% for QE, up to 100% for pion production) 
between the predictions of various generators presently 
being used by neutrino experiments Q. That is why 
the new, quite complete data on charged Q and neu- 
tral Q pion production in charged current (CC) neu- 
trino scattering pion production cross sections obtained 
by the MiniBooNE collaboration are so important and 
could be crucial in verifying our theoretical understand- 
ing not only of this process, but also of QE scattering. 

There were early theoretical attempts to describe 
neutrino-induced pion production on nuclei before the 



1 MiniBooNE tuned their generator for the calculation of the 
stuck-pion background to their data on charged and neutral pion 
production in charged current (CC) neutrino scattering. 
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data were available 181-11011. Later there were studies of 
the measured lw/QE ratios, when the cross sections for 
pion production itself were not yet available. As is un- 
derstood now all of these suffer from the fact that the 
QE cross section in the denominator of this ratio did not 
contain the 2p-2h contributions present in the data; a 
notable exception of this is the work by Martini et al. 
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who were the first to realize the importance of the 
2p-2h excitations. The few detailed comparisons of the- 
ory with the new results of the MiniBooNE experiment 
seem to indicate that a similar problem as for QE scat- 
tering exists also for neutrino-induced pion production 
12l Il3| : the data were found to be considerably above 
the calculated cross sections. This raises the question of 
the compatibility of the new neutrino-nucleus data with 
the old neutrino-nucleon data and of possible in-medium 
effects on neutrino-induced pion production. 

On the theoretical side, all (with the exception of Ref. 
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]) of the presently available detailed calculations of 
pion production rely on the so-called impulse approxima- 
tion in which the neutrino interacts with only one target- 
nucleon at a time. The interactions in medium, apart 
from Paul-blocking and Fermi-motion, are assumed to be 
the same as in free space. In such a framework it is easy to 
understand that, e.g., the total absorption cross sections 
get slightly reduced in medium in comparison with the 
elementary reaction (i.e., the cross section summed up 
over all nucleons in the target as if they were free) simply 
as a consequence of Pauli blocking. Considerably larger 
effects appear for semi-exclusive final states which are 
strongly affected by final state interactions (FSI), which 
can alter the signature of the event. Thus, the correct 
simulation of both QE and pion production events re- 
quires a model that is able to describe the elementary 
reactions as well as the final state interactions. 

It is the aim of this paper to investigate, first of all, 
the compatibility of the old pion production data on pro- 
tons or deuterium with the new data obtained on nu- 
clear targets. Second, by detailed comparisons with data 
for single-pion production on a nuclear target we study 



how far a realistic impulse approximation model of these 
processes can go in explaining the pion production data. 
Any remaining discrepancies could then possibly be at- 
tributed to primary many-particle interactions. 

In section [XT] we give a sketch of the GiBUU model. In 
section IIIII wc compare our calculations, which use two 
different parametrizations for elementary pion produc- 
tion as lower and upper bounds, with the MiniBooNE 
data. In particular, we investigate the influence of FSI 
on the absolute values and the shape of the observed 
distributions. Possible origins of remaining discrepan- 
cies will then be discussed in section IIV1 Our findings 
are summarized in sectionfV] Details on elementary pion- 
production inputs are given in appendix [AJ The influence 
of the medium modification of baryon properties on the 
various cross sections is discussed in appendix [Bl 

II. GIBUU TRANSPORT MODEL 

MiniBooNE has reported the cross section for the so- 
called CC "observed single-pion production" , that is for 
events with 1 pion of a given charge (there are data for 
7r + and 7T°) and no other pions in the final state, inde- 
pendent of the initial neutrino interactions vertex. Thus, 
to describe those events theoretically, one needs a model 
with all relevant channels included: QE scattering, A and 
higher resonance production, background 1-pion produc- 
tion, and DIS; all these reaction mechanisms are con- 
tained in GiBUU. 

In this paper we, therefore, employ the GiBUU model 
for the analysis of the data. This model has been devel- 
oped as a transport model for nucleon-, nucleus-, pion-, 
and electron induced collisions from some MeV up to tens 
of GeV. Several years ago neutrino-induced interactions 
were also implemented for the energies up to few GeV. 
Thus, we can study all kind of elementary collisions on all 
kind of nuclei within the unified framework. The code is 
written in modular FORTRAN and is available for down- 
load as open source [14{. Recently the GiBUU code was 
extended to describe also the 2p-2h excitations 15| and 
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the DIS reactions [3| for neutrino-induced reactions. Be- 
low we shortly describe the nuclear model implemented 
in the GiBUU code and used in this analysis. For more 
details on the GiBUU model see the review in Ref. [l7 

A. Elementary neutrino interactions 

A necessary input into all transport codes are the el- 
ementary reaction cross sections. Neutrino and antineu- 
trino interactions with nucleons may result in several dif- 
ferent channels. At hadronic invariant masses below the 
single-pion production threshold, W < 1.08 GeV, only 
the QE processes vn — ¥ fi~p and vp — ¥ fi + n are possible. 
They are described within the standard form factor for- 
malism as outlined in our previous publications 0, . 
We emphasize, that for the axial mass we use the world 
average value Ma = 0.999 GeV [2l| in all calculations. 

At W > 1.08 GeV the single-pion production chan- 
nel opens. The biggest contribution comes from the 
A [P 33 (1232)] resonance, at least up to invariant pion- 
nucleon masses of W « 1.5 GeV. In Refs. 
used the ANL data on proton 
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20J we had 



22j to pin down the 



details of the axial coupling to the A resonance. We had 
used the ANL data as our standard input because here 
also the absolute values of the cross section for dc/dQ 2 
were given. In this paper we now also show results ob- 
tained by using the BNL pion production cross sections 
18] as input in order to explore the consistency of both 
input data sets with the recent MiniBooNE data. Fig. Q] 
shows that the data obtained in the BNL experiment are 
generally higher than those from the ANL experiment. It 
has been argued in [23j that the two data sets are statis- 
tically consistent when taking into account their individ- 
ual errors and possible errors in the flux normalization, 
oint analysis of both data sets was also performed in 



24j . In order to obtain an estimate for the systematic 
uncertainty in the calculations we fit here both of these 
data sets separately. The curves in Fig. [TJt show the the- 
oretical description as obtained from such a fit and as 
used as input in GiBUU. In app. [5] we give the relevant 



parameters. 

With increasing W, excitations of resonances with 
higher masses, followed by their subsequent decays, give 
increasing contributions. For the prominent A resonance, 
the A — > Ntt channel saturates the one-pion width at the 
99% level. The higher resonances, however, can also pro- 
duce two and more pions as well as other mesons. In the 
GiBUU code 13 of 19 resonances having 3- and 4-star 
ratings in the Review of Particle Properties by the Par- 
ticle Data Group are accounted for. The remaining 6 are 
unaccounted for because their electromagnetic properties 
as unknown. The resonance production is described phe- 
nomenologically, as presented in . Namely, for the 

higher lying resonances the vector couplings are taken di- 
rectly over from the MAID analysis 25j so that the model 
is consistent with all results for electromagnetic interac- 
tions. The axial couplings are obtained from PCAC. 

In the same W region, non-resonant background gives 
a noticeable contribution. These processes are discussed 
in [2^. The default parameters of the model are tuned 
in such a way as to describe the ANL data on pir° and 
n7r + final states in CC neutrino scattering on neutron. 
In this paper we also show results obtained by using the 
generally higher BNL pion production cross sections as 
input. Both ANL- and BNL-tuned inputs are shown in 
Fig. [lj> together with the corresponding data. 

DIS events are included by using the PYTHIA event 
generator. The produced particles are then further prop- 
agated by the FSI treatment of GiBUU. While the DIS 
picture is applicable at high or at least moderate squared 
momentum transfer Q 2 , at low Q 2 PYTHIA treats the 
interaction within the vector meson dominance model as 
the production of vector mesons and their subsequent 
decays. 

B. Description of initial nucleus 

GiBUU describes the struck nucleus as a collection of 
off-shell nucleons. Each nucleon is bound in a mean-field 
potential, which on average describes the many-body in- 
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FIG. 1. (Color online) Single-pion production cross section on proton and neutron targets obtained in the BNL (circles; 
solid curve) and the ANL experiments [19( (diamonds; dashed curve). The curves give the lower (ANL-tuned) and upper 
(BNL-tuned) boundaries on the elementary input as used in GiBUU. 



teractions with the other nucleons. This potential is 
parametrized as a sum of a Skyrme term depending only 
on density, and a momentum-dependent contribution. 
The phase space density of nucleons is treated within a 
local Thomas-Fermi approximation. At each space point 
the nucleon momentum distribution is given by a Fermi 
sphere, whose radius in the momentum space is deter- 
mined by the local Fermi momentum, which depends on 
the nucleon density. Within this picture, contrary to the 
Fermi gas model with constant Fermi momentum (the 
global Fermi gas model), the nucleon position and mo- 
mentum are correlated. This leads to a smoother momen- 
tum distribution with somewhat more strength at lower 
momenta and to smoother nucleon spectral functions. A 
more detailed description, comparison with the global 
Fermi gas model and numerical values of the parameters 
used are given in [2^], in particular Fig. 6 there. 



All calculations reported in this paper are based on the 
impulse approximation in which the incoming lepton in- 
teracts with a single bound nucleon, with the interaction 
vertex being the same as in the case of a free nucleon. The 
vertices for the quasielastic scattering, resonance excita- 
tions and background contribution are given in Q, [2c" 
those for 2p-2h excitations in 15[ and for DIS in [l 




the publications quoted we have shown that such an ap- 
proach can give a very good description of the measured 
inclusive double-differential (with respect to the outgo- 
ing muon) cross sections. Incorporated into GiBUU it 
can make also detailed predictions for knock-out particle 
spectra and multiplicities. 

C. Pion production 



The single-pion production cross section is then given 
by (cf. 0, sec. 4.2.1) 



da uA ^' x * = I dE 



d 3 p 
(2^)3 



P(p, E)f COTI da mod P PB (r, p)F„ (q ff , r) 



(1) 



Here / CO rr is a flux correction factor / corr = (k ■ p) / (k p ); k and p denote the four-momenta of the neutrino and 
nucleon momentum, respectively. The factor Ppb describes the Pauli-blocking and P describes the hole spectral 
function, in the local Thomas-Fermi model given by 



P(p,E)=g J d 3 re[p F (r)-\p\]e(E)5^E- 

nucleus 



(2) 



5 



here Pf(t) is the local Fermi momentum given by the local Thomas Fermi model. In this spectral function all effects 
of the nucleon potential are assumed to be contained in the effective mass m* 17] which depends on location and 
momentum of the nucleon. Finally, the factor F 7r (q 7r , r) in (JlJ describes the effects of all the FSI contained in GiBUU, 
as briefly described in sec. Ill Dl 

The cross section da med in (JTJ describes the elementary pion production cross section. It can be separated into 
terms depending on the number of nucleons involved in the initial interaction: 



dcr 



mcd 



j mcd 
uo lplhl IT 



Here the first term contains the standard cross section 
for lir production on one nucleon as described in detail in 
[s, 20, 3| as briefly summarized in app. [Altogether with 



the parameters used. It is proportional to nuclear density 
oc p. The second term contains production processes that 
involve 2 or more nucleons; it is at least oc p 2 . 

The superscript 'med' in is meant to indicate that 
the spectral functions of the particles involved can change 
in medium due to binding. In all calculations we use the 
impulse approximation, i.e., only the first term in (|3|). 
with in-medium changes for the width of the A resonance 
taken from the work of Oset and Salcedo (OS) 27[. The 
separation into lp-lh and 2p-2h contributions then in- 
volves a subtlety: if we use the OS result for the A width 
in the first term, in a diagrammatical approach this al- 
ready involves terms of 2p-2h nature. In sec. IIVI we will 
discuss the possible influence of higher-order (in density) 
processes. 

All our calculations are presented for a CH 2 target. 
Charged pions can be initially produced on the proton 
as well as on the carbon nucleus. In the latter case they 
undergo FSI. Neutral pions are initially produced on car- 
bon, except for a negligible contribution from DIS scat- 
tering on proton at higher energies. 



j mcd 
uo 2p27i Itt 



(3) 



I 

four-momentum and position of a particle in the nu- 
cleus. The spectral function of a free resonance is a Breit- 
Wigner distribution with a width equal to the free width 
of the corresponding baryon resonance. In lepton- nucleus 
interactions the outgoing particles are initially produced 
inside the nucleus, and, thus, are off-shell and bound in 
the mean-field potential described above. In medium, 
on one hand, the width is lowered due to Pauli block- 
ing. Indeed, consider e.g. a nucleon resonance, decaying 
into a pion nucleon pair. This decay is Pauli-blocked if 
the nucleon momentum is below the Fermi momentum. 
On the other hand, the width is increased by the colli- 
sions inside the nucleus. For example, via the processes 
AN NN and ANN NNN the A can disappear 
without producing a pion. Secondary pion production 
is also possible, namely via the process AN — > ttNN. 
These processes, as well as elastic scattering AN — > AN, 
all contribute to the collisional width. It is essential here 
to treat the collisional broadening consistently with the 
corresponding collision terms. This means, for example, 
that whenever contribution from e.g. ANN — > NNN to 
the A width are taken into account, then also the corre- 
sponding 3-body collisions have to be taken into account. 
This is, indeed, the case in GiBUU. 



In-medium modifications of nucleons and resonances 



As outlined in [2fJ, each outgoing baryon is described 
by its spectral function, which depends on its four- 
momentum and self-energy, which in turn depends on 



A theory to calculate this medium modification was 
worked out by Oset and Salcedo (OS) These authors 
evaluated the A self-energy using a many-body expan- 
sion in terms of particle-hole and Delta-hole excitations. 
Quasi-elastic corrections, two- and three-body absorp- 
tion are included in this width. In app. [B] we discuss the 
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influence of the OS broadening of the A resonance on the 
pion yields in some detail. 

The GiBUU code can be run with all the prescrip- 
tions for medium modifications mentioned above. The 
default option in GiBUU is using the OS result for the 
A width and neglecting medium modifications for higher 
resonances. The latter give only a relatively minor con- 
tribution and using the full collisional width is computa- 
tionally expensive. We will discuss the influence of the 
various treatments of spectral functions on the final, ob- 
servable results in some more detail in sec. II VI and app.lBl 

D. Final State Interactions 

In experiments with nuclear target the outgoing 
hadrons are masked by the final state interaction: pro- 
tons can rescatter in a nucleus and produce additional 
mesons or knock-out another nucleons. Pions, that were 
originally produced through a weak excitation of A res- 
onance, can be absorbed in the nucleus or converted to 
other mesons. Thus, the inevitable presence of the FSI 
washes out the true origin of the event and makes an 
experimentally observed signal different from what one 
would expect for the scattering on a free nucleon. FSI 
can decrease the cross sections as well as significantly 
modify the shapes of the final particle spectra. In the 
present study FSI will, for example, change the signature 
of initial pion production; through FSI this process can 
lead to knock-out nucleons which mimic QE scattering. 

The restoration of the true cross section of the process 
of interest, therefore, cannot be achieved by pure experi- 
mental means and crucially relies on theoretical modeling 
implemented in the event generators. For a reliable inter- 
pretation of experiments one thus needs a model which 
provides a realistic description of both initial neutrino- 
nucleus interaction and the final state interaction of the 
produced hadrons. 

In GiBUU FSI are implemented by solving the semi- 
classical Boltzmann-Uehling-Uhlenbeck (BUU) equation. 
It describes the dynamical evolution of the phase space 



density for each particle species under the influence of 
the mean field potential, introduced in the description 
of the initial nucleus state. Equations for various parti- 
cle species are coupled through this mean field and also 
through the collision term. This term explicitly accounts 
for changes in the phase space density caused by elas- 
tic and inelastic collisions between particles. For a more 
detailed discussion of FSI see Ref. [17 1. 



E. Model validation 

The GiBUU model has extensively been tested and 
has been found to give a good description of experimen- 
tal data for very different reaction mechanisms. For the 
present investigation comparisons with pion-nucleus and 
with photon- nucleon reactions are most relevant. The 
formalism used for their description is the same as that 
just described in the preceding section. Only the neu- 
trino as incoming particle has to be replaced by either a 
pion or a photon. For both reaction types, however, the 
FSI are the same. 

Tests of pion absorption on nuclei have been passed by 
GiBUU 28j . However, these total absorption cross sec- 



tions are not very sensitive to details of the pion's interac- 
tions with the surrounding nucleus. More sensitive to the 
details of the pion interactions with the nuclear medium, 
and particularly relevant for the neutrino-induced 7r° pro- 
duction on nuclei, is the double-charge exchange reaction 
of pions with different nuclear targets. GiBUU calcula- 
tions achieve a good quantitative agreement with the ex- 
tensive data set obtained at LAMPF for incident pion 
kinetic energies of 120 - 180 MeV Q- 

Closely related to neutrino-induced reactions is the 
photoproduction of pions. On the nucleon (i.e., with- 
out final state interactions) these reactions directly test 
the vector part of the pion-production vertex. On nuclei 
this reaction tests the nuclear dynamics of pion propa- 
gation throughout the nucleus. GiBUU describes a wide 
dataset for photoproduction of neutral pions 30] on nu- 



clei for photon energies up to 0.8 GeV quite well 
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FIG. 2. (Color online) Momentum distributions of the outgoing pions for inclusive 7r° production in scattering of photons of 
energies 0.2 — 0.8 GeV off D and C nuclei. The dashed curve shows the calculations for a D target, the solid one that for C. 
The data are [jjfj; they contain coherent pion production, which is not included in the calculations; it plays a role only on the 
low-momentum side of the peak. 



As an illustration that will become relevant for the 
later discussions Fig. [2] shows the pion momentum dis- 
tributions for 7r° photoproduction off D and C nuclei. 
GiBUU calculations (which were originally done in 31 1 
and are now checked with the present version of the code) 
show a good agreement with the data. In particular, the 
shape of the distribution on deuterium (which is nearly 
equivalent to that before FSI), and on C nucleus after FSI 
are reproduced very well. The figure clearly indicates a 
significant change in the shape due to the FSI. The main 
effect is a strong absorption of pions around momenta of 
0.3 GeV due to the excitation of the A resonance and 
its subsequent pionless decay through NA — > NN or 
NNA -> NNN. 



all use the collision-broadened width as given by Oset and 
Salcedo 2jJ. Using the free width instead would yield 



significantly too large cross sections at the peak position. 

Electroproduction data for pions are very sparse, but 
the HERMES data at 28 GeV and the recent JLAB data 



at 4 - 5.8 GeV are again reproduced fairly well 321 133|. 
Also the total inelastic electron scattering data are de- 
scribed very well as shown in Ref. 2fJ. There, in Figs. 9 



and 10, it is also shown, that the inclusion of the medium- 
modified width gives a better description of the double 
differential electron scattering cross section. In particu- 
lar, the QE peak is lowered and the cross section in the 
dip region between the QE and Delta peaks is enhanced 
in agreement with the data. 



We also note here that the calculations shown in Fig. [2] Therefore, all the neutrino-induced calculations in this 
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paper are performed using the expressions for the in- 
medium width of the A resonance as given in 27| . The 
general effect of this in-medium correction is again a sig- 
nificant lowering of the pion production cross section at 
its A peak and a broadening of the latter. These changes 
are thus the same as seen in the description of the pho- 
tonuclear data. They depend only on the in-medium 
properties of the A resonance and not on its popula- 
tion mode, and should thus be also taken into account 
in all calculations of neutrino-induced pion production. 
For neutrinos these effects are discussed in more detail in 
App.E 

III. MINIBOONE CC CHARGED AND 
NEUTRAL PION PRODUCTION 

In this section we compare our model predictions with 
the MiniBooNE data on pion production. 

The MiniBooNE energy flux [34{ peaks at 0.6 GeV and 
becomes very small above 1.7 GeV. We use this flux in 
all our calculations. 

In the following discussions the data from jf], Q] are 
plotted versus reconstructed neutrino energy, while the 
theoretical curves presented here are versus real energies. 
Here we neglect this difference. For the later comparisons 
it is also essential to note that all the experimental cross 
sections for tt + production were obtained with the full 
MiniBooNE flux. Thus, for positively charged pions the 
distributions with respect to muon or pion energy do not 
depend on any energy reconstruction scheme. This is not 
so for the ir° data where a cut for the (reconstructed) neu- 
trino energies was imposed: only neutrino energies be- 
tween 0.5 and 2.0 GeV were taken into account [7]. Since 
this cut is not possible without a generator-based recon- 
struction procedure the ir° data may thus contain some 
model dependence. We also note that all our earlier cal- 
culations published in Refs. Q, [l3, 
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used the full 



MiniBooNE flux, without this cut, for all charge states. 
The net result of using now the neutrino energy-window 
for 7T° is an increase of the calculated flux-averaged cross 



sections for neutral pions compared to the earlier calcu- 
lations. 



A. Energy-dependence of cross sections 

We start our discussion with a comparison and dis- 
cussion of energy-unfolded cross sections for pion pro- 
duction. The cross sections for single-pion production 
versus neutrino energy for CC l7r + and 1tt° production 
are shown in Fig. [3l 

1. 7r + production. 

Figs. [3k,b,c show the results for lir + production. Im- 
mediately noticeable in panel (a) is that the BNL ele- 
mentary input leads to a significantly larger cross section 
than the ANL input also for a nuclear target (CH2 in this 
case), simply reflecting that fact that already the elemen- 
tary BNL cross sections are about 30% larger than those 
from the ANL experiment (see FigQJ. In the following 
we consider the cross sections obtained with the BNL 
and ANL input as upper and lower boundaries, respec- 
tively, that illustrate the systematic uncertainties in the 
final theoretical results. For both inputs the calculated 
cross section after FSI (solid lines) is about 30% lower 
than the corresponding cross section before FSI (dashed 
lines). For the cross section after FSI, which is the one to 
be compared to the data, the slope of the upper bound- 
ary (BNL input) is about the same as the slope of the 
data (within experimental errors). The lower boundary 
(ANL input), on the other hand, is flatter than the data. 
As a consequence its difference to the data increases with 
energy. 

Fig. shows the origin of the lir + events before FSI. 
Most of them come from initial A resonance production 
and its following decay (dash-dotted line). This chan- 
nel is dominant up to about 0.9 GeV incoming neutrino 
energy. Some events are background ones (dotted line). 
There is also a small contribution from higher resonances, 
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FIG. 3. (Color online) Integrated cross section for l7r + (a,b,c) and l7r°(d,e,f) CC production versus neutrino energy. Panels 
(a) and (d) contain the results of calculations using ANL- (lower curves) and the BNL-tuned (upper curves) elementary input 
both before and after FSI. Panels (b) and (e) show the contributions of various event origins to the calculated cross sections 
before FSI. Panels (b) and (f) show the same after FSI. Here the BNL-tuned input has been used. Data are from [(J, El- 



and at E v > 1.2 GeV DIS processes start to play an in- 
creasing role. 

Fig. [3J: shows the origin of the lir + events after FSI. 
As can be seen by comparing panels (b) and (c) , FSI no- 
ticeably decrease the A-originatcd single-pion production 
due to the absorption iVA — > NN; the similar process is 
possible also for other resonances. Once a pion is pro- 
duced, independent of its origin, it may also undergo a 
charge- exchange ir + n ir°p process, which depletes the 
7r + channel as the dominant one. The ensuing reduction 
of the l7r-background channel due to FSI is also notice- 
able. Other possibilities for pions to disappear include 
ttN -> ujN , (j>N, UK, AK. 

A minor amount of pions comes from the initial QE 



vertex (long-dashed line), which is only possible due to 
FSI, when the outgoing proton is rescattered. Here the 
main contribution is from the pN — > N'A — > N'N it 
reaction. Other possibilities to create pion during the 
FSI would be ujN ttN, <f>N -> ttN, ttN -> ttttN. 



When comparing to experiment, one sees that the 
lower boundary (ANL input) lies considerably (rs 40 — 
60%) below the data with the discrepancy increasing with 
neutrino energy. For the upper boundary (BNL input), 
on the other hand, the discrepancy to the data becomes 
significantly smaller. However, even this theoretical cross 
sections is still by about 15-20% below the data. 
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2. 7T° production. 

Fig. [3ji shows the results for l7r° production. Here 
the curves before and after FSI are not so different. 
This is mainly due to the fact, that side feeding from 
ir + n — > tt°p, which decreased the charged pion out- 
put, simultaneously increases the neutral pion output. 
The reverse process gives only minor relative contribu- 
tions, because the initial ir° production cross section is 
around 5 times lower. This is, however, partly compen- 
sated by charge exchange to the ir~ channel through 
7r°n — > ir~p and other channels mentioned discussed 
above for charged pion production. 

Figs. |3l3,f show the origin of the lir° events before and 
after FSI, respectively. Here the A channel is less impor- 
tant than in the case of charged pions, while background 
processes, and at higher energies higher resonances and 
DIS play a relatively larger role than for 7r + . Compar- 
ing these calculations with experiment one sees that the 
lower boundary (ANL input) lies again about 40% below 
the data at 1 GeV. This discrepancy becomes smaller 
at higher energies due to a change in slope at about 1.2 
GeV, which comes from the opening of higher resonance 
excitations and DIS events. The upper boundary (BNL 
input) also shows this bend, but is much closer to the 
data (within about 20% at 1 GeV, and within the exper- 
imental uncertainties at higher energies). 

Our predictions for 7r~ production are shown in Fig. [4] 
Without FSI the production of ir~ is negligible, a few 
events may only come from higher mass resonances in 
the processes like R + — > pp° — > pir + 7r~ or R + — > pA° — > 
ppix~ or from DIS (the curves are not shown because they 
would be indistinguishable from zero) . Cross sections for 
both l7r~ (solid curve), which is defined as one 7r~ and 
no other pions, depends significantly on the elementary 
input used; its lower to upper boundaries (ANL to BNL 
input) differ by a factor of 5. The uncertainty band for 
multi-7r _ production, which is defined as events with at 
least one it~ and any number of pions of other charges, is 
significantly narrower. Events with more than one pion 
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FIG. 4. (Color online) Integrated cross section for CC n~ 
production versus neutrino energy. The solid curves give the 
results for 1 n~ and no other pions, the dashed curves those 
for any number of pions, including at least 1 tt~ . In each case 
the upper curve gives the results obtained with the BNL input 
and the lower ones those with the ANL input cross sections. 

show up only at neutrino energies above ~ 1 GeV. The 
cross sections for tt~ production appear to be around 10 
- 20 times smaller than those for ln + . However, with 
the current statistics of the MiniBooNE experiment they 
should be measurable. 



B. Lepton observables 

The MiniBooNE collaboration has recently also pub- 
lished distributions versus muon kinetic energy, da/dT^, 
muon angle with respect to the neutrino beam direc- 
tion, da I 'd cos 8 M , and squared four-momentum transfer, 
da/dQ 2 [6[. While the former two distributions in prin- 
ciple do not depend on energy reconstruction the latter 
distribution does depend on it because Q 2 has to be re- 
constructed. There is, however, a subtle difference here 
in the way how the data for ir + and 7r° production were 
obtained. Whereas the former indeed are averaged over 
the full neutrino flux as given in 34{ the data for ir° 
were actually restricted to the incoming (reconstructed) 
neutrino energy range from 0.5 to 2.0 GeV and thus con- 
tain some model dependence. In our calculations we have 
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FIG. 5. (Color online) Kinetic energy distributions of the outgoing muons for the and liv production in the MiniBooNE. 
Data are from [6|, LZ|- The dashed curves give the results before FSI, the solid curves those with all FSI included. 
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FIG. 6. (Color online) Angular distribution of the outgoing muons for the ln + and production in the MiniBooNE. Data 
are from jj. The curves are as in Fig. [5] 



taken this restricted energy range into account by nor- 
malizing the experimental flux distribution to 1 over only 
this limited energy ranged 

In Figs. [H[l)]and [7] we compare these data with our cal- 
culations averaged over the published MiniBooNE flux. 
As in the previous subsection, each cross section is calcu- 
lated with ANL- and BNL-tuned elementary inputs, the 
corresponding curves should be considered as upper and 
lower boundaries of the theoretical prediction. For each 



2 In our earlier calculations :13j we had used the full MiniBooNE 
flux, without any restrictions. This led to even lower cross sec- 
tions. 



distribution the results are given before (dashed lines) 
and after (solid lines) FSI. 

As can be seen from these figures FSI hardly influence 
the shape of the distributions versus muon observables. 
In particular, in the cos# M distribution (Fig. [5]) there is a 
slight suppression at forward angles and an enhancement 
at backward angles; correspondingly the Q 2 (Fig. [7]) dis- 
tribution becomes somewhat flatter. The shape of the 
muon kinetic energy distributions (Fig. [5]) stays practi- 
cally the same. This insensitivity is due to the fact that 
the only effect of FSI on these muon observables is that 
they can remove events in which an initially produced 
pion (or A) was later on reabsorbed and bring in events 
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are as in Fig. \5\ 

in which the pion was produced only during FSI. 

The curves after FSI are now compared with the ex- 
perimental data. For 1tt + production, when referring 
to the shape-only comparison, our curves correspond to 
the data. The muon kinetic energy distribution is con- 
sistently 20% (upper boundary) to 60% (lower bound- 
ary) lower than the data, see Fig. [5] For Q 2 distribution 
in Fig. [7] the lower boundary is consistently 60% lower 
than the data, while the upper boundary is around 20% 
lower and a bit (within experimetal errors) flatter than 
the data. 

For l7r° production, the shape of the curves evidently 
differs from the shape of the data. In some kinematical 
regions (low T^, forward angles, moderate Q 2 ) our theo- 
retical bands are below the data, while in other regions 
(moderate T M , backward angles, low Q 2 ) they correspond 
to the data. 

The lower bounds (obtained with the ANL input) af- 
ter FSI and even before FSI are significantly lower than 
the experimental data for all three distributions for both 
lir + and lir° events. For example, for ir + the calculated 
der/dT^ cross section after FSI amounts (at the peak) 
to only about 65% of the measured value. Even before 
the sizeable pion-FSI the calculated cross section is below 
the data. For ir° events, even before FSI the same cross 
section amounts to only 65% of the data; the one after 




production in the MiniBooNE. Data are from [g, |7(. The curves 

FSI amounts to 55%. In this case the FSI have a rela- 
tively small effect. This is due to charge exchange that 
counteracts the pion absorption in this channel. In line 
with this finding is the observation that FSI only slightly 
change the shape of the curves for lir° events. 

C. Pion Observables 

Fig. [5] presents the results of our calculations for the 
kinetic energy distribution of the outgoing pions. Similar 
to the muon-related distributions, the lower boundaries 
of our cross sections are significantly lower than the data. 
The upper boundaries again gives results closer to the 
data. However, for both charge states the shape of the 
calculated and measured distributions is very different; 
this disagreement we will discuss below. 

Fig. [5] shows the pion angular distribution of the of 
the CC lir + and l7r° production. Here data are available 
for 7T° only; for forward scattering they are significantly 
higher than the lower boundary of our calcualtions, and 
only slightly higher compared to the upper boundary. 
The forward peaking is reproduced, but it is not as strong 
as exhibited by the data. 

Unlike the case of muon-related observables, for pion- 
rclatcd distributions the FSI noticeably change their 
shape. Surprisingly, the shape before FSI (see Fig. [5]) is 
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FIG. 8. Kinetic energy distribution of the outgoing ty + and momentum distribution of the outgoing it for single-pion production 
at MiniBooNE. Data are from [(j, Q]. The curves are as in Fig. [5] 
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FIG. 9. Distribution of the outgoing ty + and n° in their angle relative to neutrino beam for single-pion production at MiniBooNE. 
Data are from |7|]. The curves are as in Fig. [5] 



similar to that observed in the data. However, the shape 
of the calculated distributions after FSI looks markedly 
different. In particular, there is a significant lowering 
around T„ « 0.2 GeV for ir + and at around pi around 
0.3 GeV, as a direct consequence of the AttN dynam- 
ics in nuclei. The following processes are important for 
this structure: pion elastic scattering in the FSI decreases 
the pion energy, thus depleting spectra at higher energies 
and accumulating strength at lower energies. Simultane- 
ously, there is charge exchange scattering. At the same 
time pions are mainly absorbed via the A resonance, that 
is through irN — s- A followed by AN — > NN, which 
leads to the reduction in the region of pion kinetic energy 



0.1 — 0.3 GeV. For ir° production the additional increase 
of the cross section at lower energies comes from the side 
feeding of the 7r° channel from the dominant ir + channel 
due to the charge exchange scattering ir + n — > ir a p. In- 
verse feeding n°p — > n + n is suppressed, because at the 
energies under consideration, about 5 times less 7r°s than 
7r + s arc produced. The change of the shape of the spec- 
tra due to FSI is similar to that calculated for neutral 



current In production in 35 [ 



The particular shape calculated here for the neutrino- 
induced pions is in line with that observed experimentally 
in (7, 7T°) production on nuclear targets Q (cf. Fig. [2J. 
Since the shape depends on FSI and since the FSI are 
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the same in both neutrino-induced and photon-induced 
reactions the absence of this special shape in the neutrino 
data is surprising. 

The channel ir°n — > n~p is important for ir~ produc- 
tion in the FSI. Distributions for ir~ events are predic- 
tions of the GiBUU calculations; for completeness they 
are shown in Fig. [101 It is noticeable that the angular 
distributions for 7r~ are considerably less forward-peaked 
than for the other pion charge states. This reflects the 
fact that nearly all of the ir~ mesons are created by FSI. 

IV. DISCUSSION OF RESULTS 

Here we now summarize the comparison with experi- 
ment (cf. Fig. [3]). For the ANL input, already the results 
for the energy unfolded total 1tt + production cross sec- 
tions without FSI lie clearly (by « 20%) below the data; 
those with FSI included are about 40-60% below the ex- 
perimental data. For the BNL input, the latter discrep- 
ancy still amounts to about 15-20% for all energies, close 
to experimental uncertainties. 

For l7r° production one observes a similar discrepancy: 
the lower boundary (ANL input) after FSI is consider- 
ably lower than the data; only for the highest energies 
the calculated curve lies within the error bars. The up- 
per boundary (BNL input) is lower than the data for 
neutrino energies below 1.15 GeV. At higher energies 
the theoretical curve is within experimental uncertain- 
ties and even rises with energy a bit steeper than the 
data. 

It is interesting to note that the neutrino generator 
NUANCE, before any tuning, also gives cross sections 
that are consistently below the data. This can be seen, 
e.g. in Figs. 20 - 23 in and Figs. 8 - 14 in Q. Here for 
7r + the calculations account only for about 80% of the 
measured values, whereas for 7r° the discrepancy is much 
larger; for this channel only about 60% of the experi- 
mental value is obtained. This is qualitatively consistent 
with the results presented here. 

For the differential flux averaged cross sections as a 



function of T„ and of 9 n a more detailed picture emerges. 
Here it is seen that for both charged and neutral pion ki- 
netic energies up to about 60 MeV the data are very 
well reproduced for both charge states. However, above 
this kinetic energy the calculated cross sections are sig- 
nificantly lower than the data. As discussed before the 
lowering of the cross section in this region is due to ab- 
sorptive processes through the A resonance. It has been 
observed in photonuclear reactions and its absence in 
the neutrino-induced pion production is, therefore, un- 
expected. 

An uncertainty contributing to the remaining discrep- 
ancy could be a too low neutrino flux assumed by Mini- 
BooNE in obtaining these pion data. In 3J] the flux and 
energy spectrum have been carefully determined. Dif- 
ferent from other experiments the flux has been obtained 
from a GEANT simulation of the primary hadronic beam 
production cross sections for pions and kaons and their 
decay into neutrinos. The remaining uncertainties are 
estimated to be of the order of 10 %. In the analysis of 
the quasielastic data there is indeed some indication that 
the flux assumed by MiniBooNE may be too low by this 
amount [4] . A flux correction by this amount would bring 
the upper boundary of our calculations (BNL-tuned in- 
put) into agreement with the data for ir + , except for the 
region around « 0.2 GeV, where the special shape 
discussed above is not present in the data. For 7r° the 
calculations are with this flux renormalization within the 
experimental errors (see Fig. 1111) . 

We note that other theoretical calculations and in par- 
ticular comparisons with the MiniBooNE data are ex- 
tremely sparse. We have already pointed out that the 
results obtained with the NUANCE generator are qual- 
itatively similar to our results obtained with the ANL 
input: a large underestimate of the measured 7r° yield 
and a much smaller underestimate for tt + . We are not 
aware of any other comparisons of CC pion spectra in the 
literature. For NC there exists a comparison 37] with re- 
sults obtained both with a tuned version of NEUT and 
with a result obtained within NuWro [38f both of which 
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FIG. 10. Distribution of the outgoing n versus their kinetic energy and their angle relative to the neutrino beam for the 
MiniBooNE. The upper curves gives the result with the BNL input, the lower one that with the ANL input. 




FIG. 11. Left panel: Kinetic energy distribution of the outgoing ir + . Right panel: momentum distribution of the outgoing 7r°. 
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Data are from 6, 7], but multiplied with a factor 0.9. 



describe the data. While the 'inner workings' of the for- 
mer are not known to us, the latter model uses an un- 
conventional implementation of formation times for pion 
production through the A resonance and thus leads to 
an underestimate of pion absorption. Finally we men- 
tion the calculations of Martini et al. [ll| . These authors 



calculated in an advanced model that included RPA cor- 
relations and many-particle interactions also total pion 
cross sections. While the calculated values seem to be 
of the right magnitude a detailed comparison is not pos- 
sible since these calculations contain no final state in- 
teractions. The latter also is true for a recent approach 
extending the superscaling approximation to pion pro- 



duction _ 

Summarizing this section we conclude that the ln + 
data obtained by MiniBooNE seem to suggest that the 
higher elementary BNL data for pion production are cor- 
rect and that the ANL data underestimate the elemen- 
tary production cross section. After a possible flux renor- 
malization by about 10% the calculations based on the 
BNL data are compatible with the MiniBooNE data on 
pion production. A discrepancy with theory still exists 
regarding the pion momentum distributions. Here the 
calculations, in agreement with experimental results on 
photoproduction of pions on nuclei, predict a suppression 
in the pion spectra at about the A resonance which is not 
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clearly seen in the published neutrino data. 

A. Many-body mechanisms in pion production 

All the possible explanations discussed so far still use 
the impulse approximation in which the incoming neu- 
trino only interacts with one nucleon at a time. Many- 
body interactions are only taken into account through 
the collisional broadening of the A resonance. The latter 
is well established in photo- and electro-nuclear reactions 
and has, therefore, also been included here. If the actual 
reaction mechanism is more complicated than this and in- 
volves explicit many-body interactions then a component 
is missing in our calculations. In the experiment, the en- 
ergy reconstruction would be affected and lead to wrong 
results. A similar situation has indeed been found to exist 
for quasielastic scattering where events identified exper- 
imentally as quasielastic scattering instead cont ained a 
significant contribution of 2p-2h events [ll|, \l^L 



40|-|42| 



A corresponding mechanism leading to pion production 
would be, e.g., the process v + N + N — s- /i + A + A with 
the A finally decaying into Nir. Such processes are not 
explicitly included in our calculations. The 2p-2h con- 
tributions to QE scattering become relevant in the dip 
region and under the A resonance. It is, therefore, natu- 
ral to expect a similar process, with one outgoing nucleon 
replaced by an outgoing A to become relevant at energies 
that are about 300 MeV higher than those where 2p-2h 
becomes relevant for QE. Room for such a process could 
thus be in the region of larger energy transfers beyond 
the A excitation. It is, however, a theoretical challenge 
to cleanly separate these processes from successive scat- 
tering in a transport description in order to avoid double 
counting. 

Another process that could lead to larger pion cross 
sections is the production of off-shell pions. From de- 
tailed studies of the tt — N — A dynamics it is well known 
that the pion potential in nuclear matter is for larger 
momenta attractive. Thus, one could expect an increase 
in initial (before FSI) pion production. However, first, 



the lowering of the pion mass goes together with an in- 
crease of the A mass; thus any pion production through 
the resonance would be suppressed. Second, the off-shell 
pion has to be propagated out of the nucleus. Studies 
of the effect of pion selfenergies on the pion production 
cross sections in heavy-ion reactions have given the result 
that the final, observable pion yi eld is h ardly affected by 
the presence of a pion potential 
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453 because in dense 



matter most pions are absorbed and the last interaction 
of the pion happens for rather low densities p < 0.5po 
where the selfenergies are small. For neutrino-induced 
reactions on nuclei this point remains to be investigated. 

An estimate of an upper limit for the importance of ex- 
plicit many-body effects could be obtained from a closer 
inspection of the pion photoproduction data in Fig. m 
where the theoretical prediction has been obtained with- 
out any explicit many-body production mechanisms; only 
the collisional width of the A resonance has been used. 
This figure shows that the total cross section obtained is 
about 20% too low at photon energies around 500 MeV. 
This number thus sets a limit to any possible many- 
particle production processes, assuming that their im- 
portance is similar for vector and axial couplings. 

V. CONCLUSION 

In this paper we have calculated the neutrino-induced 
pion production on a nuclear target and have compared 
the results with the recent MiniBooNE data. The calcu- 
lations are based on the impulse approximation in which 
the incoming neutrino interacts with only one nucleon at 
a time. We have found that the experimentally measured 
total cross sections are always higher than the theoretical 
ones. However, there exists a strong dependence of the 
calculated result on the elementary cross sections used 
as input. The latter are not well constrained; the higher 
data set (from the BNL experiment) gives total theoret- 
ical cross sections close (within ss 20%) to the data. 

An inspection of the pion spectra shows that this dif- 
ference can be attributed to the pion momentum distri- 
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butions which in the theoretical calculations show a clear 
effect of pion absorption through the A resonance; this 
effect is less pronounced in the data. 

A final comparison, in which the experimental flux has 
been renormalized by 10%, gives results for the differen- 
tial cross sections that are generally in quite good agree- 
ment with the MiniBooNE data, when the BNL elemen- 
tary data are used as input. The discrepancy in the pion 
spectra at a pion momentum of about 0.3 GeV is still 
there and somewhat more pronounced for ir + production. 
However, in both cases the calculated results lie at the 
lower end of the experimental error bars so that the dis- 
agreement is statistically not well established. This near 
agreement relies on the use of the in-medium collision- 
broadened width of the A; without it the cross sections 
would be significantly larger in the peak region and the 
spectral shape would be worse. 

This leaves us with the unsatisfactory conclusion that 
there are two possibilities to explain the data. 



a remeasurement of the elementary cross sections would 
clarify the situation. With reliable elementary cross sec- 
tions at hand one would have a clear handle on the im- 
portance of many-body effects in neutrino-induced pion 
production on nuclei. 
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Appendix A: Pion production through the A 
resonance 



• First, the BNL data describe the elementary cross 
sections correctly. Then, together with a flux renor- 
malization within the experimental uncertainties, 
the data obtained by MiniBooNE can be explained 
without invoking any many-body production pro- 
cesses. Without flux renormalization many-body 
effects could amount to about 10 - 20% at most. 

• Second, The ANL data describe the elementary 
cross sections correctly. In this case the nuclear 
production cross sections are clearly too low and 
many-body production effects would have to be in- 
voked to bring the calculated cross section up to the 
measured values. A change of the flux within its 
experimental uncertainty would not be sufficient. 
Instead, many-body effects would be needed to gen- 
erate about 40% of the cross section. 

Even though the first alternative sounds more convinc- 
ing to us, in the present situation it is not possible to 
make a definite choice between these possibilities. Only 



For easier reference we collect here briefly the essentials 
of our calculation of the pion production cross section 
dv'ipihi tt m Eq. ©. For further details we refer to Ref. 

0. 

The cross section for the resonance excitation in the 
reaction £(k)N(p) — > £'(k')R(p') is in general given 
da R |k'| A{p 12 ) 
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Here the spectral function A is given by 
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The free resonance width T is momentum dependent; it is 



obtained from the Manley analysis 46J with a pole value 
of 0.118 GeV. 



The matrixelement in (jAll) is given by 
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where 9c is the Cabibbo angle. The hadronic tensor for 
resonance excitation on one nucleon, Hr, is 
1. 
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Here A is, for the A resonance, the spin 3/2 projector. 
Vertex T of weak nucleon-resonance interactions involves 
phenomenological form factors, T = 7°r^7°. Its vector 
part is completely determined by electron scattering. Its 
axial part, on the other hand, contains the crucial form 
factor generally denoted by C$ (for more details see sec. 
II. B. 2 in 20|). The resonance contribution to the pion 
production cross section dai p ih i„ is then obtained by 
multiplying daR with the free branching ratio for decay 
into ttN and assuming an isotropic angular distribution 
in the A restframe. 

The form factor C§ has to be determined by a fit to 
elementa ry pion production data. The ANL data for 
da/dQ 2 
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22J have been fitted in [20| by 
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with a = - 1.21 and b = 2 GeV 2 , with an axial mass 
of M^(ANL) = 0.95 GeV and Cf(0) = 1.17, consistent 
with PCAC. 

We have now fitted also the BNL data with the same 
functional form. The fit parameters differ from those for 
the ANL data only by the axial mass which is larger: 
M^(BNL) = 1.3 GeV. 

While these fits determine the resonance contribution 
the necessary background terms have been determined 
following again [20] . There the vector part a^ G was again 
determined by a fit to electron data, using the MAID 
analysis as a basis. For the axial part (including the 
vector-axial interference) it was then assumed, for sim- 
plicity, that it is proportional to the vector part, so that 
the total cross section then is given by 



dfTBG = (1 



(A6) 



The total pion production cross section for the \plh l7r 
final state is then given by a sum of resonance and back- 
ground contribution where the latter also includes the 
resonance-background interference term 

I\ 



A fit to the ANL data was achieved with b v * (ANL) = 3 
and b nw+ (ANL) = 1.5. The fit to the BNL data has given 
the values 6^° (BNL) = 6.0 and 6 n7r+ (BNL) = 3.0. 

All calculations in the present paper have been done 
with these parameters. The curves in Fig. [1] contain 
not only the resonance and background contributions 
described here, but also contributions from higher reso- 
nances and DIS; the latter start to contribute for neutrino 
energies above about 1 GeV. 

Appendix B: Influence of medium modification of 
the Delta 

In all of the calculations in this paper we have used 
a medium modification of the A width according to the 
Oset and Salcedo (OS) model [27j. In medium the width 
r in Eq. (|A2[) acquires a collisional contribution which is 
given by 

2 



r co ii — — - 



I p ■ 



=ImS co n(p ) , 



(Bl) 



where the selfenergy £ co ii is a function of density and 
momentum (see Eq. (4.4) in [271] )■ In the neutrino reac- 
tions investigated here the momenta of the A resonances 
produced can range up tow 1.5 GeV, i.e., well above the 
range of validity of the OS model where pa < 0.3 GeV. 
As a consequence, use of this A width requires a signif- 
icant amount of extrapolation for the higher momenta. 
For simplicity, we freeze the collisional width at the high- 
est calculated value. The price one has to pay for this 
is a loss of spectral strength. At the higher momenta 
about 10 - 15% of the spectral strength is missing (see 
Table 7.5 and discussion in 47J). All calculations in this 



oipih iff = a R-p, r obg 

1 tot 



(A7) 



paper, except otherwise noted, have been done using this 
collisional width. 

Here we illustrate in some more detail how this col- 
lisional broadening influences the cross sections. Since 
the OS model involves many-body effects on the A colli- 
sional width it is essential that the parametrization used 
is consistent with the actual collisions in GiBUU. This is 
indeed the case; in particular, when using the OS model 



19 



then also 3-body collision terms are automatically turned 
on. 

The influence of initial state interactions and medium 
modifications of the A resonance on the A production 
cross section for neutrino scattering off carbon nucleus is 
illustrated in Fig. [12] The cross section for 6 free neu- 
trons and 6 free protons, all at rest, is shown as dash- 
dotted curve labeled "6p+6n". It provides a reference 
for the following discussions. The nuclear effects, namely 
the Fermi motion, the binding nuclear potential and the 
Pauli blocking of the outgoing proton, decrease the cross 
section by around 5%, which is shown as solid curve la- 
beled " 12 C". Taking into account the OS modification 
of Delta properties leads to an additional decrease of 
the cross section by another 5 — 8%. This is shown as 
dashed curve labeled " 12 C, OS" . This decrease is a conse- 
quence of the shift of strength from the A peak position 
to larger masses due to the increased broadening. This 
extra strength at the higher masses is, however, cutoff by 
the form factors so that the net effect is a lowering of the 
cross section. 



the case, but the effect is larger: the cross section is de- 
creased by 15 — 20%. This additional suppression is due 
to the fact that the A now has an increased collision 
width so that it undergo A + N — > NN before the pion 
decay has taken place. 

While Fig. [T^] shows the results for the initial A pro- 
duction, which does not depend on FSI, we now con- 
sider the influence of FSI on the charged pions, which 
originates from the A resonance production. Fig. 1131 
shows the kinetic energy energy distribution of these tt~*~ . 
Before FSI the distribution has a rather broad peak at 
T « 0.1 GeV, which is decreased by about 5%, if the OS 
modification of the A properties is taken into account. 

After FSI, the OS modification of the A properties de- 
creases the pion distribution by about 20%. This means, 
that before having a chance to decay to a pion, a signif- 
icant part of As is absorbed in the nucleus. The overall 
effect is a decreased pion production cross section which 
increases the number of pionless events and thus serves 
as a source of the fake QE-like events. 
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FIG. 12. (Color online) Cross sections for A production, 1-7T- 
production originating from A and 07r-production originating 
from A with and without the OS medium modification [2^ . 

Since the free A primarily decays to a pion and a nu- 
cleon, one would expect a similar decrease for the single- 
pion cross section. As one can easily see, this is indeed 



FIG. 13. (Color online) Kinetic energy distribution of 7r + 
produced in neutrino scattering off carbon through the weak 
production of the A resonance and its following decay. The 
neutrino energy is E u = 1 GeV. The curves labeled OS were 
obtained using the in-medium collisional width of the A from 
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